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The properties of the "standard" type Ic supernova 19941 from 
spectral models 
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ABSTRACT 

The properties of the type Ic supernova SN 19941 are re-investigated. This object is often 
referred to as a "standard SN Ic" although it exhibited an extremely fast light curve and un- 
usually blue early-time spectra. In addition, the observations were affected by significant dust 
extinction. A series of spectral models are computed based on the explosion model C021 
(Iwamoto et al. 1994) using a Monte Carlo transport spectral synthesis code. Overall the den- 
sity structure and abundances of the explosion model are able to reproduce the photospheric 
spectra well. Reddening is estimated to be E{B — V) = 0.30 mag, a lower value than previ- 
ously proposed. A model of the nebular spectrum of SN 19941 points toward a slightly larger 
ejecta mass than that of C021 . The photospheric spectra show a large abundance of iron-group 
elements at early epochs, indicating that mixing within the ejecta must have been significant. 
We present an improved light curve model which also requires the presence of ^^Ni in the 
outer layers of the ejecta. 
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1 INTRODUCTION 

With the proposed connection between core-collapse supernovae of 
type Ic (SNe Ic) and long-duration gamma- ray bursts (GRBs; see 
[Piran 1999 for a review), interest in SNe Ic has recently increased. 
The proposed picture of SNe Ic is the collapse of the core of a mas- 
sive star which shed its hydrogen and helium envelope prior to the 
explosion. In contrast to type la supernovae, core-collapse super- 
novae exhibit a very diverse appearance. The sub-class of type Ic 
supernovae which have been occasionally connected to GRBs are 
characterised by large kinetic energies in their ejecta and are of- 
ten referred to as "hypernovae" or broad-lined supernovae. The 
progenitors of this class of objects are proposed to be very mas- 
sive stars with masses above 25A/q and are likely to form black 
holes after their collapse, while SNe Ic originati ng from lower mass 
stars p robably leave a neutron star behind fe.g.. lNomoto & Umedal 
l2002li. 
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In order to understand the properties of SNe Ic in general and 
how "normal" SNe Ic relate to the class of "hypernovae," a system- 
atic analysis of well-observed objects is needed. SN 19941 is re- 
garded as_a_grototypicalSN Ic. It has extensive observational cov- 
erage JYokoo et alJll994 frilir)Denko etalJll995t iRichmond et alJ 
1996; Clocchiatti et al. 1996), and it has been the subject of vari- 
ous theoretical studies. Both synthetic light curves a nd spectra have 
been published previou sly ( Nomot{^^lJJl994|Jlwamot^^lJ 
11994 baron et alJll996t JMillard et aLlll999t iBaron et alJll999l) : 
however, so far no consistent picture of the properties of this ob- 
ject has emerged. 

SN 19941 is also important because it seems to represent the 
branching p oint of the model sequ ence of core-collapse supernovae 
proposed bv lNomoto et al.l (2003). Therefore, it could provide the 
key for a more systematic understanding of the evolution of the 
progenitors of the various flavours of core-collapse supernovae. 

Ilwamoto et al.iri994) used a low-mass model (C021) to fit the 
light curve of SN 19941. Their model has only ~ 1A/q of CO-rich 
ejecta, suggesting that SN 19941 was the stripped C-0 core with a 
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mass of ^ 2.1Mq originating from a ~4Mq He-core. This corre- 
sponds to a star witli a main sequence mass of ~ ISAf© . To achieve 
such a thorough stripping, repeated interaction with a binary com- 
panion must be envisaged. 

In this work we present a series of spectral models for the 
early-ti me spectra o f SN 19941 base d on the explosion model 
C021 by llwamotoe t al. ( 1994); see also lNomoto et alJ<1994 . This 
explosion model provided an excellent fit to the light curves of 
SN 19941. However, light curve models alone are not able to con- 
strain the explosion parameters uniquely because the light curve 
shape depends on a combination of kinetic energy and ejecta mass. 
The peak luminosity is constrained by the amount of ^''Ni synthe- 
sised during the explosion. To determine the total luminosity, it is 
necessary to have reliable values for parameters such as reddening 
and distance. In addition, the bolometric correction used to infer the 
luminosity has to be derived from the shape of the spectral energy 
distribution. Using the highly parametrized s pectral code SYNOW 
JFisheretai]ll999l. and references therein), iMillard et all <1999h 
provided line identifications for some of the spectra of SN 19941. 
[Baron et al. ( 1999) have calculated non-local-thermodynamic equi- 
librium (non-LTE) models for a few spectra based on the outcome 
of C021 and concluded that this model is not able to reproduce 
the observed spectra in detail. They suggested that a higher mass 
model would be required to match the line features better. 



2 SPECTRAL MODELS 



We used the observed spectra published by Filiooenko et al. 
tl99^) which have bee n flux calibrated against the photometry of 
iRichmond et alJ il99a) . Our series of photospheric models cov- 
ers the observations from April 4 to May 4, 1994. (UT dates are 
used throughout this paper.) For all models the explosion date was 
set to be March 27, 1994 (JD 2 44 9 439), correspo nding to a rise 
time to B-band maximum of 12 d jlwamoto et alJll994') . In what 
follows the epochs are given in days relative to this date. The B- 
band maximum of SN 19941 occurred on April 7 (JD 2449 450) 
t Rich mond et a l. 1996). We h ave excluded the first spectrum (April 
2) in Filippenk o et alj il995h because the wavelength coverage of 
this observation is too narrow for a conclusive model. Spectra later 
than May 5 already show too much nebular emission to be treated 
with a photospheric approach. In addition, we used a spectrum from 
April 3 that was taken at the Multiple Mirror Telescope (MMT) by 
Schmidt et al. (1994, priv. comm.). 

The early-time spectral models have been computed employ- 
ing the method which w as first propos ed by Abbott & Lucv ( 1985) 
and further developed bv lMazzali & Lucv> 1 1993), Lucv ( 1999), and 
IMazzalL L200Q') . The method employs a Monte Carlo simulation of 
the line transfer based on the Sobolev approximation and includes 
a line-branching scheme. At the end of the Monte Carlo calculation 
the emergent spectrum is derived from a formal integral solution 
on the basis of the source fu nctions t hat can be extracted from the 
Monte Carlo simulation (see lLuc\il999.') . 

For the model computation the ejected material is divided into 
an outer region, where the true continuum opacity is assumed to be 
negligible, and an inner region which is not treated explicitly. At the 
dividing boundary of the computational area a thermal blackbody 
continuum is imposed with a temperature that is adjusted according 
to the luminosity of the emergent radiation at the outer boundary. 
Note that this temperature has to be determined in an iterative pro- 
cess because it is a priory not clear what fraction of the radiation 
is scattered back into the core region. This constitutes the effect of 



line blanketing where the back-scattered radiation causes the mat- 
ter to heat up. The ejecta are treated in spherical symmetry. At this 
stage we assume a homogenised composition for the entire compu- 
tational volume. 

To fit the early- time spec tra we used the density structure from 
the C021 explosion model ( Nomot o et al J [19941 llwamoto et alj 
[l994) . The model is expanded homologously to match the desired 
epoch. The composition was adjusted in order to obtain a good fit to 
the observed spectrum. The abundance distribution from the nucle- 
osynthetic yields of the C021 model averaged over velocity space 
above each epoch's photosphere served as a guideline. We adopt a 
distance to the host galaxy of 8.4 Mpc, correspon ding to a distance 
modulus of /i — 29.6 mag JFeldmeier et alll997h . 



2.1 Reddening 

The amount of extinction by dust has a significant impact on the 
estimate of the total luminosity of the object and, therefore, the 
amount of synthesised ^^Ni derived from models. T he extinction 
toward SN 19941 has b een a controversial topic ( Baronet alll996t 
IRichmond et alJll99a) and could not be firmly established. The 
presence of a strong, narrow Na I D line in the spectra, as well as 
the colours, suggest that the object was significantly reddened, but 
the exact amount is subject to large uncertainty. While the analysis 
of the Na I D line i mplie s a very large extincti on of -E(-B — V ^) ~ 1 
liHo & FiliDpenkcl ll995l) . spectral models bv iBaron et alj il99d) 
suggest a value for E{B — V) between 0.30 and 0.45 mag. 

The amount of reddening affects the model spectra because 
the luminosity affects the temperature and hence the properties of 
the spectrum (given a fixed distance and density structure). The 
temperature structure determines the overall shape of the spec- 
trum and, via the ionisation balance, the strengths of characteristic 
line features. A test series of different values for E{B — V) be- 
tween 0.15 and 0.50 mag showed a preferred range between 0.25 
and 0.35 mag. For larger reddening the ionisation tends to be too 
high to reproduce th e characteristic lines of singly ionised species 
( Millard et al. 1999). A value of E(B~ V) near the lower end re- 
sults in temperatures that shift the line ratios, in particular those of 
Ca II, too much toward the transitions with low excitation energies. 
For our full series we chose E{B — V) = 0.30 mag, however note 
that the uncertainty for this value may be as large as ±0.05 mag. 
For models with larger reddening the composition required to fit 
a spectrum contains more iron-group elements at the expense of 
lighter elements because a larger reddening requires a larger tem- 
perature, shifting the spectrum to the blue. This blue-shift has to be 
counterbalanced by stronger line blocking, which is efficiently pro- 
vided only by iron-group elements having a large number of lines. 

In the following the properties of each model are briefly dis- 
cussed. 



2.2 Model series for the photospheric epoch 

2.2.1 t = 7d, Aprils, 1994 

For this spectrum the calibration against the corresponding pho- 
tometry required a correction of 0.17 mag to match the R and V 
bands reasonably well. However, with this correction the U and B 
bands are still too bright by 0.67 and 0.30 mag, respectively. This 
deviation is too large to apply a single averaged value. Therefore 
we re-calibrated the spectrum to match only the V and R bands, 
which represent the majority of the flux. Keeping this inconsistency 
in mind, we did not try to match the observations exactly but rather 
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Figure 1. Synthetic spectra of SN 19941 in comparison with the observed spectra. The observations are from lFilippenko et aljil99j) . except for the spectrum 
of April 3, which was taken at the MMT (Schmidt et al. 1994, priv. comm.). For each spectrum the epoch in days after the assumed explosion date (March 27, 
1994) is indicated. The flux F\ is given in units of 10"^'* ergcm"^ s~^ A^. In panel c the ions that form the major contribution to the observed features are 
indicated. The identifications do not change significantly at different epochs. 



adjusted the model to obtain the basic features and reach agree- 
ment with the U , B, V and R photometric points for this epoch. 
Consequently, the fit does not look too good at all wavelengths (see 
Fig. ^), but the main features are reproduced, except for the O I 
A7774 line, which is significantly weaker in the observed spectrum. 
This line tends to be too strong in all our models. It is a highly sat- 



urated line that is fairly insensitive to changes in the oxygen abun- 
dance, within reasonable limits. 

Apart from the inconsistency in the data, part of the problems 
to fit the overall slope of the observed spectrum could also arise 
from the model assumption of a thermal lower boundary which may 
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not sufficiently represent the physical conditions in the ejecta. We 
discuss this possibility in Sec. 13. 21 



2.2.2 t^8d, April 4, 1994 

The total luminosity required for this model is i — 
1.89 X 10*^ergs~^. With a photospheric velocity of Vp^ = 
11 600 km s~^ this results in a temperature of the underlying black- 
body of Tbb ~ 9679 K. In the blue and visible wavelength bands 
the fit is reasonable, while the slope of the pseudo-continuum of 
the model does not match that observed in the redder parts of the 
spectrum. The spectrum shows various strong absorption features 
which can be attributed to the ions indicated in Fig^. Interest- 
ingly, the absorption at ~ 6100 A is only partially produced by 
Sill; it also contains a significant contribution from various Nel 
lines, most important Ne I 6402. As before, the absorption feature 
of the O I A7774 line appears too deep. 

The model also does not reproduce the spectral slope of the 
observation redward of about 6500 A. Part of this problem may, as 
in the previous spectrum, be due to a calibration problem in the 
data. The correction that was applied to match the V and R pho- 
tometry results in an offset of 0.14 mag in the I band. Therefore, 
the /-band magnitude of the model is actually in closer agreement 
with the photometry than the integrated magnitude of the spectrum 
in that bandpass filter. Nevertheless, also for this case the possibil- 
ity of a physical effect that is not correctly represented in the model 
should be kept in mind (see Sec. l3.2t . 

Another interesting feature of this spect rum is the absorption 
structure between 4000 A and 4200 A which lMillard et alj <I999h 
attribute to Sc II. In our models we can reproduce that feature using 
significant amounts of Sc, but this also causes a variety of other 
Sc lines to show up elsewhere in the spectrum. In the later spectra 
the observed absorption seems to grow into a distinct feature which 
is mostly caused by Till and Fell lines, while there is no clear 
indication for Sc lines anymore. Therefore we suggest that even in 
the early spectrum the main components of this feature are Fe lines 
that are not correctly reproduced by our model. In addition, other 
ions such as C II and Si II contribute to this feature to some extent. 



2.2.3 t = 10 d, April 6, 1994 

This spectrum (Fig.^) represents the obse rvation closest to the B - 
band maximum, which occurred on April 7 ( Richmond et al. 1996). 
In this model, the continuum slope is already better reproduced than 
in the first spectra although there is still an offset redward of ~ 
7000 A. We used a photospheric velocity of Wph = 10800 kms"^. 
The total luminosity increased to L — 2.31 x 10^^ ergs^^, and 
the temperature of the incident radiation was Tbb ~ 9230 K. With 
the exception of the "Sc" feature and red continuum slope, the fit is 
reasonably good. The observation starts to show the separation of 
the Fe II feature at 4900 A which is reproduced by the model at this 
epoch. Also, the depth of the Ca II infrared lines is not sufficient in 
the model while the Ca II H&K absorption is fitted well. 



2.2.4 t = IQ d, April 12, 1994 

This spectrum (Fig.Ql), 16 d after explosion (corresponding to 5d 
after B-maximum), exhibits clear differences compared to the pre- 
maximum spectra, and it is significantly redder. The parameters that 
have been used for this fit are L = 2.10 x lO''^ erg s~^ for the total 
luminosity and a photospheric velocity of Uph — 8800 km s ~ ^ . The 



resulting temperature at the inner boundary was Tbb = 7690 K. 
The double structure between 4500 and 5000 A is more pro- 
nounced, and our model mostly fails to reproduce it. The feature 
is mainly formed by Fe II lines that are all part of the same multi- 
plet (Fell AA4923, 5018, 5169). The fact that we see the lines sep- 
arated in the observation indicates that either the density or at least 
the ionisation fraction of Fe II must have a steep gradient in the ve- 
locity region just above the photosphere, which is not reproduced 
by our models of the earlier epochs. At later epochs the models 
resolve this structure better. 

Apart from that, most of the features in the model are sys- 
tematically blue-shifted compared to the observation. Even with 
a drastic change of parameters this problem remained, suggesting 
that the density structure is too flat within the velocity range where 
most lines form at this epoch. 



2.2.5 t = I'd d, April 15, 1994 



The luminosity used for this model (see Fig. ^) decreased fur- 
ther to 1/ = 1.46 X 10^^ erg s^^ at a photosphere located at 
iiph = 7200 km s"^ with a temperature of Tbb = 6999 K. In this 
spectrum the double feature of Fe lines below 5000 A also starts to 
show in the synthetic spectrum. 

Starting from this epoch the slope of the near-IR flux is well 
reproduced by the model. This may be surprising because at later 
epochs the approximation of a photosphere becomes less accurate. 
Given the compactness of model C02I, it may be that the photo- 
spheric approximation holds reasonably well even at epochs well 
after maximum (see i|3.2t . 

The model does not provide clear evidence for the origin 
of the deep absorption observed at 10 500 A. In previous work 
iFiliDDenko et al."l995VClocchiatti et al.''l996VMillard et al."l999t 
iBaron et al. 1999), various suggestions for the origin of this feature 
have been discussed, but no conclusive answer could be reached. 
We will investigate this feature in il3.3l 



2.2.6 t ^ 2Q d April 16, 1994 

This spectrum (Fig. ^) is just a day after the previous one and 
therefore fairly similar. A model with L = 1.30 x 10'*^ergs~^, 
Wph = 6800 km s~^, and Tbb = 6774 K reproduces most features 
quite well, except that some of the peaks tend to be somewhat low. 
The narrow peak at 6560 A, which is not reproduced by our 
model, could to some extent be due to emission of Ha by a sur- 
rounding shell. To test this possibility, however, is beyond the scope 
of this work. 



2.2.7 t = 2Ad April 20, 1 994 

The i = 24 d epoch represents the last of the early spectral series. 
Overall this spectrum is still similar to the one taken 4 days earlier, 
except that the Fe II and Na I lines become deeper. For the model in 
Fig.^ we used L = 9.18 x lO''^ ergs~^ and a photospheric ve- 
locity of Uph = 5400 kms~^. The temperature of the photosphere 
was derived to be Tbb = 6383 K. 



2.2.8 t^38dMay4,1994 

This spectrum, 38 days after explosion, is already quite late to be 
treated with a photospheric code because it may contain a signifi- 
cant fraction of nebular-line emission mixed into the emission com- 
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Table 1. Nebular models for different values of reddening. 



4000 5000 6000 7000 8000 9000 10000 
rest wavelength [A] 

Figure 2. Nebular model for the nebular spectrum 159 d after explo- 
sion compared to the observation assuming a reddening of E{B — V) = 
0.30 mag. F\ is in units of 10~^^ ergcm"^ s~^ A^ 



ponents of the P-Cygni profiles. Therefore, it is not surprising that 
some of the emission features cannot be reproduced by this ap- 
proach (e.g., the Call IR triplet). Nevertheless, the fit is reasonable 
(see Fig.0i). The luminosity used was L = 3.49 x 10'*^ ergs~^; 
the photosphere was set to «ph = 2100 km s~^. These values are 
not strongly constrained because a clear definition of a photosphere 
is not possible at such late epochs. Rather, the velocity sets the ra- 
dius which is required to obtain the correct temperature structure to 
reproduce the correct ionisation ratios given the chosen luminosity. 
The temperature at the photosphere was Tbb = 6129 K. 



2.3 Nebular Spectrum, t = 159 d 

As is typical for SNe Ic, the nebular phase develops early because 
of the small ejecta mass. To test the inner part of the ejecta we also 
modelled the nebul ar spectrum taken on September 2 using a non- 
LTE nebular code ( Mazzali et alj |200ll) which assumes a constant 
density, one-zone model inside a given velocity. 

The spectrum has an epoch of 159 d after the proposed ex- 
plosion date. It was c orrected based on the photometric data of 
[Richmond et alj il99dl but in the blue part the correction is some- 
what uncertain because of the presence of an underlying galaxy 
continuum. However, this uncertainty does not significantly affect 
the conclusions in this section. 

For the model shown in Fig.|2|we assumed an expansion ve- 
locity of Wexp = 5500 km s^^ at an epoch of i — 159 d. The 
model has a mass of M = O.43A/0 below this velocity. To re- 
produce the correct line strengths using the same distance and red- 
dening values as for the photospheric spectral series, 0.07A/q of 
^^Ni are required. This generates a luminosity of the nebula of 
L = 2.22 X 10"° ergs"^ with a faction of ~ 10 per cent of the 
7-rays depositing their energy in the gas. The remaining composi- 
tion needed to reproduce the spectrum includes O.22M0 oxygen, 
0.09 Mq carbon, O.O2M0 silicon, and O.O2M0 sulphur. 

To test the impact of reddening we calculated a series of mod- 
els with different values for E{B-V): 0.15, 0.30, and 0.45 mag. 
Within this range of E{B — V) it is possible to find a set of param- 
eters that fit the observed spectrum well, but the masses of the vari- 
ous models are different. The resulting parameters for these models 
are summarised in TableQ 

The mass enclosed within 5500 kms^^ is larger than the cor- 
responding mass of C021, which is only 0.1 A/0. Also, the line 
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Figure 3. Mass fraction of the various elements versus the velocity of the 
photosphere ("Uph) that was used for the respective model. At each epoch 
the models are calculated assuming a homogeneous distribution of the ele- 
ments. 



profiles do not suggest the presence of a mass cut at a velocity of 
~ 2000 kms^^ but indicate the presence of material at lower ve- 
locities. 



3 DISCUSSION 

3.1 Composition and derived magnitudes 

Fig.|3|shows a plot of the mass fraction for some elements used for 
the respective models versus the photospheric velocity, Uph. Note 
that we assume a homogeneous composition for each model; hence, 
this plot will only roughly represent a slice of the ejecta compo- 
sition. In particular, for features that form significantly above the 
photosphere, the velocity may not be representative. Also, for el- 
ements that do not contribute any distinct line features, the abun- 
dance is not strongly constrained. 

An interesting aspect of this plot is that the abundance of iron- 
group elements (Ti, Cr, Fe, Co, and Ni) decreases toward lower 
photospheric velocities. This is a feature which is seen in hyper- 
novae that exhibit strong jet-like asymmetries. The fact that we 
find a higher concentration of heavier elements in the outer part 
of SN 19941 may suggest that the explosion of this object also 
showed asymmetry to some extent, although the nebular spectra 
do not show a strong indication for that. 

In comparison to the original C021 model, we get a signif- 
icantly larger mass fraction of oxygen at low velocities at the ex- 
pense of iron-group elements. The high abundance of iron-group 
elements at velocities above ~ 7000 km s~^ is also not predicted by 
C021 even if we assume some contribution of primordial material 
with solar composition added into the outer part of the explosion 
model. The abundance of Si is significantly lower than predicted 
by C021, while C and Ne seem to be more abundant. On the basis 
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Figure 4. Comparison of observed (squares) and model (crosses) magni- 
tudes in different wavelengtli bands. The observed magnitudes for B, V, 
and R have been obtained by interpolating the photometry points from 
iRichmond e t al. ( 1996) to the observation dates of the spectra. The fZ-band 
points are from Schmidt, et al. (1994, priv. comm.). 



of the spectral models it is not possible to put strong constraints 
on the detailed composition structure; in particular, for elements 
that do not have many spectral lines, the error bars may be quite 
large. Nevertheless, there is strong indication that the composition 
is more strongly mixed than in the theoretical ejecta model. 

Fig.|4|shows the magnitudes derived from the models in com- 
parison to the "observed" magnitudes. The data points for B, V, 
and R have been obtai ned by interpolating the original points of 
IRichmond et alj i 19961) to the epochs of the corresponding spec- 
tral observations. The L'^-band points are taken from Schmidt, et 
al. (1994, priv. comm.) and are shown at the original epochs be- 
cause there is insufficient spectral coverage of this wavelength 
band. Since the models resemble the observed spectral shapes rea- 
sonably well, the magnitudes from the models are overall in good 
agreement with the observed ones. Only at the later epochs the U- 
band magnitudes are underestimated by the models. 



3.2 The location of the "photosphere" 

Fig. |5| shows the photospheric velocities (tiph) derived from the 
models as a function of epoch. Also indicated are the measure- 
ments of the line velocity from the observed spectra assuming that 
the A6100 feature is due to Si II 6347, 6372. Our models suggest 
that there may be a significant contribution from Ne I 6402. This 
would result in an underestimate of the "real" velocity because the 
Ne I line has a slightly redder rest wavelength. The upper and lower 
limits correspond to different assumptions on the blend of lines that 
makes this feature. The lower dotted curve corresponds to a fit to 
the feature with a single Gaussian profile, while the upper dotted 
curve represents the estimate from a fit with two Gaussian profiles 
assuming that the bluer one is Sill. At the early epochs the pho- 
tospheric velocities of the models clearly exhibit a different trend. 
Therefore, the absolute velocities deduced in this way may not be 
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Figure 5. Photospheric velocities of the models as a function of epoch in 
comparison to the velocities measured from the "Si 11" feature at 6100 A 
(see text). 



a good representation of the photospheric velocities and should not 
be over-interpreted. 

Because the photospheric velocity of the models are primarily 
constrained by the temperature that is needed to match the shape of 
the spectrum and the degree of ionisation, this different behaviour 
of the Si II velocities may result from the failure of the approxima- 
tion to describe the shape of the early-time spectra with a thermal 
blackbody continuum. The mismatch between the observed and 
model red continua give another indication that this description is 
not appropriate for those epochs. Generally, in hydrogen-deficient 
supemovae the assumption of a thermal continuum photosphere is 
questionable because the physical conditions present in the ejecta 
do not result in sufficient true continuum opacity to achieve ther- 
malization (e.g., Sauer 2005, Sauer et al, in prep.). Imposing this 
thermal continuum in the models generally leads to an offset of the 
model flux in the red and infrared wavelength bands. 

This problem is less severe in SNe Ic than in SNe la be- 
cause the composition by mass fraction is dominated by the lighter 
species (oxygen and neon) while in SNe la the iron-group elements 
make the largest contribution. This means that even for the same 
mass density, the number density in SNe Ic is larger by a factor of 
up to 2. Since the true processes (bound-free and free-free transi- 
tions) that are needed to thermalize photons scale with the number 
density of the ions and free electrons, the processes are more ef- 
fective in SNe Ic. In addition, the stronger influence of Thomson 
scattering (which is also proportional to n^) will result in a larger 
number of scattering events for each photon before it can escape 
the ejecta. This effectively increases the volume in which the pho- 
tons have a chance of being thermalized. This may explain why the 
spectral shape of SNlc spectra can be described fairly well with a 
thermal continuum to later times whereas in SNe la, the approxima- 
tion to assume a blackbody continuum at the photosphere becomes 
increasingly worse for later times. 

The early-time spectra of SN 19941 are very blue suggesting 
that they may have a stronger component of non-thermalized radi- 
ation than the later ones. This may be due to a steep gradient in the 
density in the outer parts of the ejecta. If a substantial quantity of 
7-ray photons is generated in the outer, thinner layers of the ejecta, 
the effective path lengths of the photons may become too short to 
allow the photons to be thermalized. A prerequisite for the presence 
of a non-thermalized radiation field at early times is, however, a sig- 
nificant amount of radioactive ^''Ni located in the outermost layers 
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of the ejecta. This outer layer of Ni would not strongly contribute to 
the later light curve, as a large fraction of the 7-ray photons would 
escape without depositing their energy in the ejecta. Some ^ Ni is 
indeed needed in the outer part of the ejecta to reproduce the early 
light curve rise (see S|4} 

For a more consistent investigation of this phenomenon a 
better description that includes the non-thermal formation of the 
pseudo-continuum is needed. 

3.3 The infrared feature at 10 300 A 

Spectra that extend into the near infrared (Fig.^-h) show a dis- 
tinct and broad absorption feature aroun d 10 300 A. Originally , 



this feature was attributed to He I AlO 830 ipilippenk o et al. 19^ 



IClocchiatti et al. 1996) . However, as pointed out by Baron ej 



tall 



11999t) . this identification would require a detection of at least He I 
A5876 at similar velocities (~ 15 000kms~^) in the observed 
spectrum. To model the He I lines correctly it is necessary to con- 
sider the non-thermal excitations by fast electrons caused by Comp- 
ton scattering of 7-photons from the ^''Co decay. Without those 
processes the excitation of He I at the temperatures present in t he 
ejecta is too low to cause any significant absorption iLucvll99in . 

Since our code does not treat non-thermal processes we can- 
not address this problem in a self-consistent way. To nevertheless 
test the possibility that He I causes the 10 300 A absorption we 
enhanced the optical depth in the He I lines introducing by hand 
a "non-thermal factor" that mimics the effects of departure from 
LTE caused by the fast electrons produced by the deposition of 
7-photons fsee lHarkness et alll987l : lTominaga et aL , 2005). We al- 
lowed this factor to vary with velocity to obtain a good fit of the 
10 300 km s~^ absorption. Although some lines arise from the sin- 
glet system and others from the triplet system, such as A6678, we 
a pplied a singl e en hancement factor to all He I lines. The studies 
of lLucvl <199lh and lMazzali & Luc\J (1998) find that the departure 
from LTE is comparable for both states. We chose the observed 
spectrum of i = 19 d (cf. Fig.^) for this study because that is the 
earliest observation that extends far enough to the IR to show this 
feature. In the first model shown in Fig.|6^ we enhanced the optical 
depth in all He I by a factor that varied from 10* at 16 000 km s~^ 
to 10^^ in the outermost shells above 30 000kms^^. The abun- 
dance of helium was uniformly set to 10.5 per cent by mass. Most 
of the strong He I lines that also show up in the model spectrum 
blend in with other features in the spectrum, however the model 
inevitably shows a clear absorption at ^ 5500 A. AlO 833 is the 
strongest He I line in our model and thus may show up before any 
other lines can be clearly identified. Unfortunately, the observed 
spectra that cover the infrared wavelengths have a gap in the region 
where the He I A5876 lines would be expected. The observed spec- 
tra of the later epochs that do cover this wavelength range show a 
small absorption dip in the peak around 5450 A, which could be 
identified with the He I A5876 lines, but the feature is not clearly 
present in the early epoch spectra. Apart from this, the observed ab- 
sorption is much weaker than what the model spectrum would sug- 
gest for this line. Also, the corresponding velocity for He I would 
be 26 000 km s~^, which is inconsistent with the identification of 
He! AlO 833 for the AlO 300 absorption. Of course, the method 
adds more free parameters and therefore the abundance of He I can- 
not be constrained but depends on the choice of the enhancement 
factor (for a given density of the explosion model). The very large 
enhancement factors needed to obtain a reasonable fit, however, in- 
dicate that an unreasonably large amount of several solar masses of 
helium would have to be present at the high velocities to obtain this 
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Figure 6. Models for the t = 19 d spectrum. In the upper panel (a) 
the strength of the helium lines was ailificially enhanced above v Ri 
16 000 km s~^ by factors between 10** and 10^^ to mimic non-thermal 
excitation in an outer helium-rich shell (assuming a homogeneous mass 
fraction of 10.5 per cent for helium). It is not possible to reproduce the 
broad absorption feature at 10 300 A in this way without also increasing 
the strength of the He I A5876 line to a level that is not observed. The labels 
in Panel a indicate the positions of the strongest He I lines at a velocity of 
vhc = 17000 kms~. The middle panel (b) shows a model obtained by 
increasing the optical depth in the C I line by factors of 10^ to 2 X 10'^ above 
14 000 km s~ ^ to fit the 10 300 A absorption. The mass fraction of carbon 
was 21.7 per cent. This model shows a variety of C I lines in other places 
that are not observed. The positions of the C I lines con'espond to a velocity 
of vq = 9000 km s~ ^ . The lower panel (c) shows a model with both He I 
and C I enhanced. The enhancement factor for He I above 16 600 km s~^ 
reaches from lO'' to 10^^, for C I we applied a factor of 50 in the entire 
ejecta. This leads to an acceptable fit of not only the IR-absorption but also 
various other Hnes at '^ 9000 A, which are not reproduced by the original 
model (Fig.^). The flux is given in units of 10~^* erg cm~-^ s~^ A^. 
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absorption considering that to account for the non-thermal excita- 
tion alone should not require a factor of more than a few 10'' tLucy, 
ll99lU1Vla77,a1i&T,uc^^l99S^ . 



Millard et alJil993) argue from their highly parametrized line 
fits with S YNOW that the 10 300 A feature could originate from a 
blend of He I, Si I, and C I lines. In our models of these epochs the 
ionisation equilibrium is generally such that Si II is more abundant 
than Si I by several orders of magnitude. Therefore, the mass in Si 
needed to obtain considerable absorption in Si I is very high and 
leads to a significant overestimate of the Si II lines in the spectra. 
The strongest line of Si I in the observed wavelength range is Si I 
AlO 790. If we assume that this line is primarily responsible for the 
observed IR absorption, it follows that the velocity of this Si I layer 
should be around 15500 kms~^. In that case one would also ex- 
pect the presence of absorption features by the much stronger Si II 
AA6241, 6349 lines but none of the spectra suggest the presence 
of such a high-velocity silicon layer. Given those indications, we 
regard the identification with Si I as unlikely. 

The presence of carbon at velocities below 11 000 km s~^ at 
day 19 is consistent with the earlier spectra at i = 8 d and 10 d: The 
small features in the peak around 4100 A in those early spectra (see 
line label in Fig.Q;) can be partially attributed to ClI A4267 at a 
velocity around 11 000 km s~^, although our models do not repro- 
duce the details of the spectrum in that region The smaller absorp- 
tion features between the Ca II IR-triplet and the deep absorption at 
10 300 A in the later spectra can be identified with C I lines. These 
lines are indicated in Fig.lSJ), which shows a model in which the 
optical depth in C I was enhanced above 14 000 kms~^ by factors 
between 10^ and 2 x 10^ to fit the 10 300 A absorption. The mass 
fraction of carbon in this model is 21.7 per cent. The positions of 
those lines correspond to a velocity of around 9500 kms~^, which 
also agrees with the velocity of the O I lines near 7500 A. The en- 
hancement factors needed for C I are significantly lower than for 
helium. In addition, the observed spectra of earlier epochs are con- 
sistent with the presence of C II lines at similar velocities. Never- 
theless, the model shown in Fig.|6j5 does also not provide a reason- 
able fit to other spectral regions. That suggests that carbon alone 
may also not be the cause of the 10 300 A absorption. 

The lower panel in Fig. |6| shows a model with different fac- 
tors mimicking non- thermal excitations in both He I and Cl. This 
combination gives a good fit not only to the IR-absorption but also 
to the C I-lines shown in Fig. |6j). In addition, the amount of He I 
needed is low enough not to cause strong absorptions elsewhere 
in the spectrum. For this model the optical depth in He I was en- 
hanced by factors between 10^ to lO'^^ above 16 600kms~^. For 
C I a uniform factor of 50 was applied in the entire ejecta. The rel- 
atively small enhancement-factors needed for C I and the fact that 
one would not expect that non-thermal effects are significant for 
this ion, suggest that a slightly different composition and ionization 
structure at larger velocities may be responsible for the formation 
of the C I lines. 

For a more conclusive determination the effect of varying 
composition as well as a more physical description of the non- 
ther mal excitations h as to be considered. Also, as pointed out by 
e.g. IPrzvbilla i2005l) . in full non-LTE calculations of early type 
stars the strength of the IR-He I line sensitively depends on the cor- 
rect treatment of line blocking as well as the details of the chosen 
model atom. Therefore, fully consistent non-LTE models will be 
necessary to model the formation of this line. 



3.4 Oxygen abundance and total mass 

Common to all synthetic fits of the photospheric epochs presented 
in SJ2|s the fact that the predicted O I A7774 line that is too strong. 
This line is highly saturated and does not depend sensitively on the 
abundance assumed for oxygen. Therefore, one has to adopt a very 
low oxygen abundance to fit the observed depth of the O I line. In 
the t = 19d spectrum, for example, the oxygen mass fraction has 
to be reduced to ~ 4 per cent. At the same time, however, a signifi- 
cant oxygen mass of 0.22A/q is needed in the innermost region to 
reproduce the nebular spectrum (see i|2.3t . 

Another possibility is that the total oxygen mass fraction is 
indeed larger, but that, owing to non-thermal ionisation processes 
the ionisation of O I is actually higher than derived by our models. 
O II does not have many lines from low-lying levels in the observed 
wavelength bands. The presence of unexplained absorption dips in 

the spectra around 7000 A could hint at forbidden O II lines, which 

I 1 1 I 

have also been suggested by Millard et al. ( 1999). Non-thermal ion- 
isation can be induced by some amount of mixing of ''"Ni-rich ma- 
terial into the oxygen layers. The larger amount of oxygen in the 
inner region seen in the nebular phase supports this scenario. 

The obvious contradiction in this picture, however, lies in the 
identification of C I lines in the infrared as discussed in the previous 
section. If C I is seen at velocities similar to that of the O I line, 
it is not clear why non-thermal ionisation processes should affect 
oxygen but not carbon. 



3.5 Bolometric light curve 

The bolometric light curve is one of the common ways to constrain 
a theoretical model. Ideally, to obtain a bolometric light curve from 
observed data one uses a set of well-calibrated spectra from UV 
through IR wavelengths. The available data for SN 19941 does not 
cover the JHK bands in the infrared. We computed the bolomet- 
ric light curve of SN 1 9941 using the photometric data in UBVRI 
JRichmond et alJ 19961 and B. Schmidt, priv. comm.). This UBVRI 
integrated light curve from SN 19941 data (open squares in Fig.0 
alone represents a lower limit to the bolometric light curve since the 
IR contribution is not taken into account. In addition, we derived an 
upper limit estimate by adding the IR-contribution of SN 2002ap 
( Yoshii et al. 2003), which is the most similar object for which in- 
frared data have been published. For this addition we assume that 
the relative contribution of the infrared bands in SN 19941 was the 
same as for SN 2002ap after the data points of SN 2002ap have 
been rescaled to the appropriate extinction. The upper limits are 
shown as crosses in Fig. For SN 19941 additionally one point 
in each of J and H and two in the K band are available from 
iRudv et al] (l99J). This allows us to compute a single bolometric 
point for the epoch f = 28 d after explosion which is indicated by 
the filled triangle in Fig.0 This single point agrees nicely with the 
upper limit light curve which justifies the choice of the SN 2002ap 
data to derive the upper limit light curve for SN 19941. It also shows 
that the IR contribution to the bolometric flux is not negligible. 
Generally, for a SN Ic the contribution of the IR-wavelength bands 
ranges from about 20 per cent before B maximum, to roughly 40 
per cent near K maxim um. In the late ph ase the IR-bands con- 
tribute about 30 per cent JYoshii et al. 2003t lPatat et all200lh . 

To compute the bolometric light curves, we converted the 
magnit udes into fluxes using the standard UB VRI Johnson-Cousin 
system (Buser & Kurucz 1978; Bessell 1983). Using V epochs for 
reference, we then interpolated to these epochs the flux in other 
bands. The monochromatic fluxes were used as low-resolution 
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Figure 7. Bolometric light curves for SN 19941. The squares denote the 
bolometric light curve obtained from UBVRI observations of SN 19941 
without correction for the unobserved JHK bands. The filled triangle in- 
dicates the point for which JHK is also known. The crosses are the data 
points obtained by correcting the SN 19941 points for the unobserved JHK 
bands by assuming the same contribution as in SN 2002ap. The results of 
the spectral models are shown as diamonds. 
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Figure 8. Synthetic light curve for the modified C021 ejecta model (solid 
line) in compaiison with the light curve obtained from the original C021 
model (dashed line). The filled circles show the bolometric light curve from 
the observations (see 83.51 . 



spectra and transformed to luminosities adopting a distance modu- 
lus of m — M = 29.6 mag, an extinctio n of i?(i? — 'i/) = 0. 3 mag 
(see i|2.H . and the extinction law of Sav age & Mathij il97St) . 

The luminosities used in our spectral models are shown as di- 
amonds in Fig.Q They are in good overall agreement with the esti- 
mate from the observations corrected for the infrared contribution. 



4 LIGHT CURVE MODELS 

Fig. m shows the synthetic light curve of the ejecta model that we 
built, as suggested by the nebular-spectrum model, by adding a low- 
velocity dense core to the density structure of C021 (solid line), 
compared with the light curve of the original C021 (dashed line) 
and with the bolometric light curve of SN 19941 (circles). We used 
the one-dimensional supernova light curve synthesis code devel- 
oped by Iwamo to et alj J2000|). The code solves the energy and 
momentum equations of the radiation plus gas in t he co-moving 
frame a nd the gray ^-ray t ransf er simultaneousl y; see lMazzali et alj 
J2002h . lDeng et alj JlOO^ . and'Mazzali' EoO^ for our approxima- 
tions on optical line opacities and Eddington factors. The ''''Ni mass 
and its distribution were adjusted to optimise the fit. For the relative 
abundances of other elements, we used the results of our spectral 
modelling. 

Our ejecta model has a total mass of ~ 1.2Af0, of which 
^OAMq is distributed below 5500 kms"^ as required in the neb- 
ular spectrum modelling (see § 2.3). For comparison, C021 has 
in total ~ O.QMq with only ~ O.IMq below 5500kms"\ The 
total kinetic energy, ~ 1 x 10^^ erg, is barely changed by the ad- 
dition of the low-velocity dense core. The best light curve fitting 
^''Ni mass is ~ 0.07Mq for both ejecta models, although the dis- 
tribution of ^®Ni is different. For C021, we filled its low-velocity 
region with as much ^®Ni as possible without spoiling the fit to the 
peak to increase the energy deposition rate at late times. This re- 
sults in ~ 0.04-0.05A/o of ^^Ni in the inner ~ O.IMq ejecta at 
V < 5500 kms~^. For our ejecta model, which has a core as mas- 
sive as ~ O.4Af0, a ^^Ni abundance of ~ 10 per cent in that core 
is enough to reproduce the slow decline of the late-time light curve 
of SN 19941, all of which is distributed above 3,000kms"^ to 
avoid overestimating the flux around day ~ 40. Both models need a 
substantial amount of '^''Ni between 8000 and 10 000 km s"^ to re- 
produce the fast-rising early light curve (~ O.OOSA/q for C021 and 
~ O.OO3M0 for our model). This is consistent with the blue nature 
of the early spectra, which indicates substantial 7-ray heating and 
^''Ni mixing in the outer part of the ejecta (see i|3.2> . 

The light curve of our ejecta model reproduces the observed 
light curve of SN 19941 well, both in the peak phase (before day 
~ 30) and in the tail phase (after day ~ 50). By contrast, the light 
curve of C021 declines too fast to follow the observations after 
day ~ 50 and the discrepancy increases to ~ 50 per cent in flux by 
day ~ 100. This confirms the results of our nebular spectrum mod- 
elling that, compared to C021, more low- velocity mass is needed 
to generate the observed line luminosity in the day 159 spectrum. 
The low-velocity dense core may suggest that the explosion of 
SN 19941 was actually aspherical. Spherical explosion models can 
not produce this dense inner core because they exhibit a cutoff ve- 
locity below which material falls back to form the compact object 
in the center of the explosion. In contrast, models of aspherical ex- 
plosions generally show a higher density in the inner, low- velocity 



core because the fall-back region is distorted (e.g. , Khokhlo^^l 
19991 iMaeda et alJl200l l2003l : iMazzali et alj |2004l: iTomita et al 



200d) . Our model light curve overestimates the flux by < 20 per 
cent at the transition phase between day ~ 30 and day ~ 50. Fur- 
ther tuning the ^^Ni distribution and the ejecta density structure in 
the low-velocity core could solve this discrepancy, which is, how- 
ever, not significant. On the other hand, this may also be a multi- 
dimensional effect, the modelling of which is then beyond the scope 
of this paper. 
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5 CONCLUSIONS 

We have re-investigated the properties of the Type Ic supernova 
SN 19941 because this object is among the best observed "normal" 
SNe Ic and possibly represents a link between different branches 
of core-collapse supernovae. The study was based on the density 
structure of the explosion model C021 by Iwamo to et alJ J 1994) 
although we allowed the composition to deviate from the predic- 
tion of this model. Overall the explosion model seems to give a 
reasonable representation in the outer parts of the ejecta, although 
the mismatch of the velocities of some features in the early phase 
may suggest that a density structure with a steeper gradient in the 
intermediate-velocity range could be more suitable. The nebular- 
phase models, in contrast, require a significantly larger mass at 
low velocities than predicted by the model. Our estimate for the 
ejecta mass inside a radius corresponding to ?; = 5500 km s~^ is 
at least O.4Af0. From the photospheric model at f = 24 d, which 
has a photospheric velocity of 5400 kms^^, one can estimate the 
mass above the outer velocity of the nebular spectrum to be about 
O.GMq. This adds up to a value that is larger than the total ejected 
mass of O.9M0 of model C021. 

The spectral models also give some indication that the explo- 
sion of SN 19941 exhibited some extent of asymmetry or at least 
inhomogeneity which was able to mix ^^Ni out into the oxygen- 
rich layers. This is also supported by the trend that the abun- 
dance of Fe-group elements actually seems to decrease toward 
the centre of the ejecta. In turn, oxygen-rich material appears to 
be present at very low velocities as can be seen in the nebular 
spectrum. In comparison model C021 does not show much mix- 
ing. Recent hydrodynamical simulations (Burrows & Haves 199a: 
iLai & Goldreich"200d':'Scheck et al. 2004) and spectropolarimetric 
observations iWang et al-2001-2003.:.Leonard et al-200Z.200a) of 
core-collapse supernova explosions suggest the presence of global 
asymmetries which could induce large mixing of the ejecta. The 
overall picture for SN 19941 is that it probably originated from the 
explosion of a low-mass core that stripped its envelope possibly by 
binary interaction. The peculiar shape of the early-time spectra, as 
well as the very fast decline of the light curve show that this object 
was by no means a "typical" SN Ic. This should be kept in mind 
when comparisons to other objects are made. 

The amount of extinction of SN 19941 by dust is subject to 
large uncertainties. Our models suggest that the value of E{B — 
V) = 0.45 mag used by Iwamoto et al. ( 1994) to model the light 
curves is probably significantly too large. A value of at most 
0.30 mag seems more suitable to fit the observed spectra. 

Given the constraints on the ejecta mass and reddening from 
the spectral models, we recalculated the theoretical light curve for 
the C021 model with a more massive low-velocity core region. 
This new model still fits the observed bolometric light curve that 
we constructed from the spectra well (including a correction for 
the unobserved infrared bands). To obtain this fit it was necessary to 
place a substantial amount of ^^Ni at velocities between 8000 and 
10 000 km s^^. This supports the corresponding conclusion from 
the spectral fits. The total amount of ^®Ni was ~ O.OTA/q and did 
not change compared to the original C02 1 model. (The effect of the 
lower reddening than assumed by Iwamot o et alJil994h is balanced 
by a larger distance.) 

The nature of the IR-absorption around 10 300 A cannot be 
revealed unambiguously on the basis of the methods chosen here. 
The most likely identification seems to be Cl at intermediate veloc- 
ities with a contribution of He I at high velocities. Neither ion seems 
to be able to account for the entire absorption without causing un- 



seen features in other parts of the spectrum. C I can also account 
for other line features which are not reproduced by our models oth- 
erwise. A more in depth analysis using a consistent non-LTE treat- 
ment of all important processes including non-thermal excitation of 
He I is needed to provide a more conclusive picture. 



ACKNOWLEDGEMENTS 

We acknowledge B. Schmidt for providing us with some data and 
S. Taubenberger for re-calibrating the spectra. This work was sup- 
ported in part by the European Union's Human Potential Pro- 
gramme "Gamma-Ray Bursts: An Enigma and a Tool," under con- 
tract HPRN-CT-2002-00294. D.N.S. thanks W. Hillebrandt at the 
Max Planck Institute for Astrophysics for his hospitality. Some of 
the authors thank the Kavli Institute for Theoretical Physics at the 
University of California, Santa Barbara for its hospitality during 
the program on "The Supernova Gamma-Ray Burst Connection," 
where this paper was completed. This research was supported in 
part by the National Science Foundation under Grants PHY99- 
07949 and AST-0307894. 



REFERENCES 

Abbott D. C, Lucy L. B., 1985, ApJ, 288, 679 

Baron E., Branch D., Hauschildt P. H., Filippenko A. V., Kirshner 

R. P, 1999, ApJ, 527, 739 
Baron E., Hauschildt P, Mezzacappa A., 1996, MNRAS, 278, 763 
Bessell M. S., 1983, PASP, 95, 480 

Burrows A., Hayes J., 1996, Physical Review Letters, 76, 352 
Buser R., Kurucz R. L., 1978, A&A, 70, 555 
Clocchiatti A., Wheeler J. C, Brotherton M. S., Cochran A. L., 

Wills D., Barker E. S., Turatto M., 1996, ApJ, 462, 462 
Deng J., Tominaga N., Mazzali P. A., Maeda K., Nomoto K., 

2005, ApJ, 624, 898 
Feldmeier J. J., CiarduUo R., Jacoby G. H., 1997, ApJ, 479, 231 
Filippenko A. V. et al. , 1995, ApJ, 450, LI 1 
Fisher A., Branch D., Hatano K., Baron E., 1999, MNRAS, 304, 

67 
Harkness R. P et al., 1987, ApJ, 317, 355 
Ho L. C, Filippenko A. V., 1995, ApJ, 444, 165 [Erratum: 463, 

818 (1996).] 
Iwamoto K. et al., 2000, ApJ, 534, 660 
Iwamoto K., Nomoto K., Hoflich P., Yamaoka H., Kumagai S., 

Shigeyama T, 1994, ApJ, 437, LI 15 
Khokhlov A. M., Hoflich P A., Oran E. S., Wheeler J. C, Wang 

L., Chtchelkanova A. Y., 1999, ApJ, 524, L107 
Lai D., Goldreich P, 2000, ApJ, 535, 402 
Leonard D. C. et al., 2006, Nat, 440, 505 
Leonard D. C. et al., 2002, PASP, 1 14, 35 
Lucy L. B., 1991, ApJ, 383, 308 
Lucy L. B., 1999, Ac&A, 345, 211 
Maeda K., Mazzali P. A., Deng J., Nomoto K., Yoshii Y, Tomita 

H., Kobayashi Y, 2003, ApJ, 593, 931 
Maeda K., Nakamura T, Nomoto K., Mazzali P. A., Patat F., 

Hachisu I., 2002, ApJ, 565, 405 
Mazzali P. A. et al., 2006, ApJ in press (astro-ph/0603516) 
Mazzali P A., 2000, Ac&A, 363, 705 
Mazzali P. A., Deng J., Maeda K., Nomoto K., Filippenko A. V., 

Matheson T, 2004, ApJ, 614, 858 
Mazzali P A. et al., 2002, ApJ, 572, L61 



© 0000 RAS, MNRAS 000, 000-000 



The properties ofSN 19941 from spectral models 1 1 



Mazzali P. A., Lucy L. B., 1993, A&A, 279, 447 

Mazzali P. A., Lucy L. B., 1998, MNRAS, 295, 428 

Mazzali R, Nomoto K., Patat R, Maeda K., 2001, ApJ, 559, 1047 

Millard J. et al., 1999, ApJ, 527, 746 

Nomoto K., Maeda K., Umeda H., Ohkubo R, Deng J., Mazzali 

P., 2003, in Hucht K. v. d., Herrero A., Rsteban C, eds, Proc. 

lAU Symp. 212, Astron. Soc. Pac, San Francisco, p. 395 
Nomoto K., Umeda H., 2002, ASP Conf. Ser. Vol 253, p. 221 
Nomoto K., Yamaoka H., Pols O. R., van den Heuvel E. P. J., 

Iwamoto K., Kumagai S., Shigeyama R, 1994, Nat, 371, 227 
Patat R et al., 2001, ApJ, 555, 900 
Piran R, 1999, Phys. Rep., 314, 575 
Przybilla N., 2005, A&A, 443, 293 
Richmond M. W. et al. 1996, AJ, 1 1 1, 327 
Rudy R., Dotan Y., Woodward C, Cole J., Hodge T., 1994, 

lAUCirc, 5991,3 
Sauer D. N., 2005, PhD thesis, Technische Universitat Miinchen 
Savage B. D., Mathis J. S., 1979, ARA&A, 17, 73 
Scheck L., Plewa T., Janka H.-T., Kifonidis K., MuUer E., 2004, 

Phys. Rev. Lett., 92, 11103 
Tominaga N. et al., 2005, ApJ, 633, L97 
Tomita H. et al., 2006, ApJ in press (astro-ph/0601465) 
Wang L., Baade D., Hoflich R, Wheeler J. C, 2003, ApJ, 592, 457 
Wang L., Howell D., Hoflich R, Wheeler J., 2001, ApJ, 550, 1030 
Yokoo R, Arimoto J., Matsumoto K., Takahashi A., Sadakane K., 

1994, PASJ,46, L191 
Yoshii Y. et al., 2003, ApJ, 592, 467 



© 0000 RAS, MNRAS 000, 000-000 



